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Abstract 
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1 Introduction 



The quark masses depend on the renormalization scheme and, within a given one, on the 
renormalization scale. The dependence on the latter is usually referred to as "running" 
and is governed by the quark mass anomalous dimension, j m . 

The two-loop anomalous dimension is known for long [1]; the three-loop result is also 
available from Refs. [2,3]. 

The quark mass anomalous dimension 7 m is defined as 

/^T^L^b = mimics) = -mJ2lial + \ (1) 



where a s = a s /7r = g 2 /(47r 2 ), g is the strong coupling constant. To calculate 7 m one needs 
to find the so-called quark mass renormalization constant, Z m , which is defined as the 
ratio of the bare and renormalized quark masses, viz. 

1 >3 



Within the MS scheme [4] the coefficients (Z m )^ are just numbers, with D = 4 — 2e 
standing for the space-time dimension. 

In practice Z m is usually computed from the vector and scalar parts of the quark self- 
energy Ey (p 2 ) and T*s(p 2 ) ■ m our convention, the bare quark propagator is proportional to 

+ Yiy{p 2 )) — mj (1 — £^(jo 2 )) . Requiring the finiteness of the renormalized quark 
propagator and keeping only massless and linear in m q terms, one arrives at the following 
recursive equations to find Z m 

Z m Z 2 = 1+K t [z m Z 2 Y? s (p 2 )} , Z 2 = 1 - K e {Z 2 Y? V { P 2 )} , (3) 



where K e f(e) stands for the singular part of the Laurent expansion of /(e) in e near e = 
and Z 2 is the quark wave function renormalization constant. Eqs. (3) express Z m through 
massless propagator-type (that is dependent on one external momentum only) Feynman 
integrals (FI), termed as p-integrals below. 

In this letter we first present the results of the analytic calculation of the quark mass 
anomalous dimension to 0{af). Then we apply it to evaluate the running of a heavy 
quark mass in QCD with the number of active flavours rif varied from 3 to 6. 
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2 Results and discussion 



Taken literally, Eqs. (3) require calculation of a host of one-, two-, three-, and four-loop 
p-integrals to find Z m to 0{af). At present a direct analytical calculation of a p-integral 
with the loop number not exceeding three is a rather easy business. First, there exists an 
elaborated algorithm — the method of integration by parts of Ref. [5,6] — which allows 
one to analytically evaluate divergent as well as finite parts of any three-loop p-integral. 
Second, the algorithm has been neatly and reliably implemented in the language FORM 
[7] as the package named MINCER in Ref. [8]. 

The situation with four- loop diagrams is quite different. At present there is simply no 
way to directly compute the divergent part of a four-loop p-integral. An indirect and 
rather involved approach to perform such calculations analytically is to use the method 
of Infrared Rearrangement (IRR) of Ref. [9] enforced by a special technique of dealing 
with infrared divergences — the -Reoperation of Refs. [10,11]. Recently, the technique has 
been dramatically simplified by explicitly resolving all the relevant combinatorics for a 
particular case of 4-loop diagrams contributing to the correlator of two scalar or vector 
currents in Refs. [12] and [13] respectively. We have extended these improvements for the 
problem at hand. 

The four-loop diagrams contributing to Eqs. (3) to order a 4 s (about 6000) have been 
generated with the program QGRAF [14], then globally rearranged to a product of some 
three-loop p-integrals with a trivial (essentially one-loop) massive Feynman integral and, 
finally, computed with the program MINCER. All calculations have been done in MS- 
scheme. As is well-known, every anomalous dimension is the same for MS- and MS schemes 
[15], thus, the results are also valid for the latter one. 

Our result for the quark mass anomalous dimension reads (in evaluating the colour factors 
we have assumed the case of the colour group SU C (N); rif stands for the number of quark 
flavours) 



7o 



N 2 ~l 
AN 



7i 



N 2 



1QN 2 



3 203 Ar2 

- H N 2 +n f 

8 24 1 



- 5 -N 
6 . 



(4) 



72 



+ n f 



N 2 - 1 I 


129 


64N 3 1 


16" 




1177 


[>- 




108 



129 Ar2 11413 , 

N H N 

16 216 



A^ 3 -6A^C(3)-6A^ 3 C(3) 



35 ,t2 

— N 

54 



(5) 



N 2 



73 



1261 50047 



256N 4 128 



+ 



384 



. , 66577, . 460151 - 
-N 2 - —^rN + N 6 + 21 C(3) 



1152 



1152 



3 



+ n f 



+ n 3 



, & + l^iV3 - ^ - ™N CIS) - 85^ C(3 ) - 8 ^ C 



l N + ±& - - ™N C (3) - 85^ 3 C(3 ) - ^ C (3) 

6 216 81 2 sv ; sv ; 6 SV ; 

+ 33iV 3 C(4) + 33iV 5 C(4) - 30N C(5) + 5(W 3 C(5) + 80iV 5 C(5)" 

1Q 8QQ 

~^N 2 + ||iV 4 + 10iV 2 ((3) + 10iV 4 ((3) - 6N 2 C(4) - 6iV 4 ((4) 

_^_iV 3 + -iV 3 C(3)' 
162 9 v . 



(6) 



Here C is the Riemann zet a- function (C(3) = 1.202056903 . . ., C(4) = tt 4 /90 and £(5) = 
1.036927755 . . .). With N = 3 one gets the following result for QCD 
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Note that in three-loop order we exactly reproduce the known result of Ref. [2,3]. The 
four- loop term proportional to n 3 is in agreement to the one found in Ref. [16]. 



In numerical form 7 m reads 



7m = - a s - a 2 s (4.20833 - 0.138889^) 

- a 3 (19.5156 - 2.28412n/ - 0.0270062nJ) 

- a 4 (98.9434 - 19.1075n/ + 0.276163nJ + 0.00579322^). 

It is amusing to compare the exact result <-y| xact = 44.2629. . . for rif = 3 and a recent 
estimation of Ref. [17] 

7 f ( nf = 3) = 7 ^ n/ = 3 } = 40.6843. (10) 
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For the case of QED with rif identical fermions our result reads 

QED 
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The RG equation (1) is usually solved by 

m(fj) c(a 3 (fj)) 
m(/i ) c(a s (/i ))' 



(14) 



where 



c(x) = (a;)' S6 {l + (7 1 -A7o)a; 

+ g [(^ ~~ &7o) 2 + 72 + AVo - /5i7i - An 

+ 7(71 - Ato) 3 + 7K71 -/3i7o)(72 + A 2 7o - A71 - /9 2 7o) (15) 
.6 2 

+\ ( 73 - /?i 3 7o + - /3 3 7b + P1I1 ~ P2I1 ~ P1I2)} * 3 + O(o: 4 )} . 



Here 7, = Ji/Po, Pi = Pi/Po, (i=l,2,3) and $ are the coefficient of the QCD beta-function 
defined as: 



i+2 



(16) 



Equation (15) explicitly demonstrates that the knowledge of the four-loop coefficients 73 
and Pz is absolutely necessary for the self-consistent running of the quark mass in the 
cases when the mass-dependent terms of order a\ are taken into account. Examples can 
be found in Refs. [12,17,18]. 



The four-loop beta-function has been recently analytically computed in Ref. [19] with the 
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result 
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Now we are in position to evaluate the c-function numerically with four-loop accuracy. 
For strange, charm, bottom and top it reads: 



Cs (x) = x 4/9 {l + 0.895062a; + 1.37143a; 2 + 1.95168a; 3 , (n f = 3), 

c c ( x ) =x 12/25 (l + 1.01413a; + 1.38921a; 2 + 1.09054a; 3 ), (n f = 4), (18) 

c b (x) =a; 12/23 (l + 1.17549a; + 1.50071a; 2 + 0.172478a; 3 ), (n f = 5), 

ct(x) =x 4/7 (l + 1.39796a; + 1.79348x 2 - 0.683433a; 3 ), (n f = 6). 
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